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Abstract 



We have found an exact formula expressing a general correlation function con- 
taining both products and ratios of characteristic polynomials of random Hermitian 
matrices. The answer is given in the form of a determinant. An essential difference 
from the previously studied correlation functions (of products only) is the appear- 
ance of non-polynomial functions along with the orthogonal polynomials. These 
non-polynomial functions are the Cauchy transforms of the orthogonal polynomials. 
The result is valid for arbitrary ensemble of /3 = 2 symmetry class and generalizes 
recent asymptotic formulas obtained for GUE and its chiral counterpart by different 
methods. 

1 Introduction 

It is a classical result of analysis known from the nineteenth century that orthogonal 
polynomials can be represented as multi- variable integrals. For example, let vrj(x) de- 
note a monic orthogonal polynomial of degree with respect to the measure dn{x) = 
exp{— NV {x))dx on the real axis. 



and TTj{x) = x^ + lower degrees. Then for the monic orthogonal polynomial tinIx) there 
exists an integral representation (see Szego [1]) 





JV 




(2) 
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In the formula above A{xi, . . . ,xn) stands for the Vandermonde determinant, and the 



N-l 



constant Zpf is given by the product Zpf — N\ H 



The classical integral representation of the monic orthogonal polynomial 7rN{x) given 
in (2) suggests a random matrix interpretation. Indeed if we consider an ensemble of 
random N x N Hermitian matrices H with the joint probability density of eigenvalues 
(see Mehta [2]) 



P(^) {x,,...,xn) = ^ exp I -TV E ^i^i) I A'(a:i , . . . , xn) 



(3) 



the monic orthogonal polynomial 71n{x) can be understood as an average of the charac- 
teristic polynomial Zj^[x, H] = det{x — H) over the ensemble 



(4) 



A natural question which arises at this point is the following: what do we obtain if instead 
of one characteristic polynomial we average a combination of characteristic polynomials 

(for example, their product). This question was recently addressed in Brezin and Hikami 
[3], Mchta and Normand [4] who derived a generalization of formula (4). The authors 
considered the ensemble average of products of characteristic polynomials. It was found 
that this average is expressed as a determinant: 



H 



A(A) 



det 



7rAr(Ai) 7rjv+i(Ai) ... TlN+L-li.^l) 
T^n{.^2) 7rAr+i(A2) ... 7riv+L-l(A2) 



(5) 



where A(A) = A(Ai, . . . , Al). (See also a related paper by Forrester and Witte [5] which 
deals with a particular case of positive integer moments: Ai = . . . = Al). 

In the present paper we obtain a further generalization of formulas (4) and (5). Namely 
we compute the general spectral correlation function of characteristic polynomials defined 
as 

L 

I n ZN\}ij,H]. 

}CM{i,il)^r4 ) (6) 

\ n ZN[ej,H]' 
i=i H 

Here the symbols e, fi denote the vectors with the components: 

e= (ei,...,eM); /i = (//i, . . . , //l) (7) 
For the function /Cjv(e, fi) to be well defined we assume that Im ^ 0, and also assume 
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N > M, N > L. Then the result we obtain for the correlation function /C7v(e, A) is 



/iAf-M(ei) hisr-M+ii^i) 



hN+L-l{^l) 



/CAr(e, p,) 




det 



(8) 



where hk{e) denotes the Cauchy transform of the monic orthogonal polynomial 7rfe(a;), 



The formula above is of interest since it shows the place of the correlation function /CAr(e, fi) 
in the theory of orthogonal polynomials. Furthermore it provides new insight on )Cn{^, A), 
for example it relates /Civ(e, /*) with Riemann-Hilbert problem for orthogonal polynomials 
(see section 5) 

Correlation functions /CAr(e, /i) are interesting objects themselves since they contain 
a very detailed information about spectra of random matrices. In particular, knowledge 
of such correlations enables to extract the n-point correlation function of spectral den- 
sities. Moreover, distributions of some interesting physical quantities such as, e.g. level 
curvatures are expressed in terms of 1Cn{^', A), (see e.g [6] and Appendix A of [7] for more 
examples) . 

The investigation of correlation functions of characteristic polynomials and of their 
moments is also motivated by hope to relate statistics of zeroes of the Riemann zeta func- 
tion to that of eigenvalues of large random matrices [8, 3] . The spectral determinants are 
also relevant for interesting combinatorial problems, sec [9]. Other numerous applications 
of spectral determinants are in the theory of quantum chaotic and disordered systems, 
and in quantum chromodynamics, see [10]-[11] for an extensive list of references. 

There are several analytical techniques for deahng with the integer moments (positive 
or negative) of characteristic polynomials and more general correlations functions. Their 
applicability varies with the nature of the underlying random matrix ensemble. The par- 
ticular case of the Gaussian measure - Gaussian Unitary Ensemble (GUE), as well as 
its chiral counterpart, chGUE - can be studied very efficiently by a modification of the 
standard supersymmetry technique [12] suggested recently by the authors [10, 11, 13]. 
One starts, as usual, by representing each of the characteristic polynomials as a Gaussian 
integral over ant i- commuting (Grassmann) variables and inverse characteristic polynomial 
by Gaussian integral over complex variables. This allows to average the resulting expres- 
sions straightforwardly. At the next step one employs the so-called Hubbard-Stratonovich 
transformation for "fermionic" degrees of freedom, and then exploits the Itzykson-Zuber- 
Harish-Chandra integrals together with their natural non-compact extensions. A detailed 




(9) 



and 



27ri 



(10) 



7n-l = 



'n— 1 
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account of the method and results as well as related references can be found in the recent 

papers [10]- [13]. 

The resulting expressions in [11, 13] revealed a very attractive determinantal structure 
which was especially evident in the case of chiral GUE (Laguerre ensemble) [13], see also 
[14]. When one deals with the positive moments only such structures naturally arises 
in the framework of the orthogonal polynomial method [3, 4]. This fact is suggestive of 
the idea that the determinantal structure in the general case (i.e. when both products 
and ratios of characteristic polynomials are involved) should be valid for an arbitrary 
unitary-invariant potential. In present paper we show that this is indeed the case. 

The paper is organized as follows. In section 2 we derive an algebraic identity which 

represents H [e^, H] as a sum over permutations. This identity together with formula 

(5) will permit us to rewrite /C7v(e, A) as an integral over M variables with subsequent 
derivation of formula (8). Section 3 illustrates our approach on the simplest case M — 1 
and L — 0. Here we derive a multi- variable integral representation for Cauchy transforms 
of monic orthogonal polynomials. The computation of the general correlation function 
/Cjv(e, /i) for any integers < L, M < is given in section 4. In section 5 we establish 
relation with the Riemann-Hilbert problem for orthogonal polynomials proposed by Fokas, 
Its and Kitaev [15, 16] (see also [17, 18]), and outhne a way to investigate the large N 
asymptotic of /CAr(e, /t) by Deift-Zhou steepest descent/stationary phase method [18, 19, 
20, 21, 22, 23]. A detailed asymptotic analysis of /C7v(e, A) is beyond the scope of this 
paper and will be given in the subsequent publication [24]. Finally in section 6 we briefly 
discuss some "duality relations" for matrix integrals emerging for the case of the Gaussian 
potential. 



2 Inverse of products of characteristic polynomials 
as a sum over permutations 

Let M < N, and xi,...,xn denote the eigenvalues of a matrix H (dim H — N x N). 
Then the following algebraic identity holds: 



M N-M 



n 



'I 



I T-r ^N-M \ 



E 

<t6 Sjv/S;,;_mXSm 



M 



1=1 Zn^.H] 

A(a:(j(i), ■ • ■ ,a^CT(M))A(a;(^(M+i), ■ • ■ ,a^(7(Ar)) 



(11) 



M 



n e, 



\i,j=l 



X 



aii) 



A(xo-(i), . . . ,Xo.(Ar)) 



where Sat is the permutation group of the indices 1, . . . , A?"; Sm is the permutation group 
of the flrst M indices and Sjv-m is the permutation group of the remaining N — M indices. 
Identity (11) follows as a consequence of the Cauchy-Littlewood formula [25] 



M N 



n 11(1 ^ = E^^(^i'---'^^)^A(yi,---,yM) 

3=1 i=l A 



(12) 
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and the Jacobi-Trudi identity [25]: 



det (^x^ 



Xj-j+N 



Sx{xi,...,Xn) = r (13) 

A(a;i, ...,Xn) 

where the Schur polynomial Sx{xi, . . . ,xn) corresponds to the partition A and the indices 
i,j take the values from 1 to A^. In order to prove (11) we rewrite the Schur polynomial 
sx{xi, . . .,xn) in (12) as 



sx{xi,...,xn) ^^^^^ 7-7 ^ (14) 

A(xi, ...,xjv) 

(As sx{yi, . . . , um) = for any partition with the number of rows larger then M only the 
partitions with Xm+i — — Xn — — ^ contribute in equation (12)). In equation (14) 
i^TT = 1 (or = 0) if the permutation tt G Sjv is odd (even). 

Any permutation vr G Sat in the above sum can be decomposed as a product of three 
subsequent permutations, i.e. n = a ■ 1:2 ■ tti- The first one, tti G Sm is a permutation of 
elements of the set [1, 2, . . . , M]; the second permutation, 7r2 G Sn-m is that of elements 
of the set [M + 1, M + 2, . . . , A?"]; and the third permutation, a G Sn/Sm x Sn-m, is an 
exchange of elements between these two sets. For example, let us take M — 3, N — 5, and 

( 1 2 3 4 5 \ 
1 4 2 5 3 /' "^^^^ permutation can be represented as follows: 

TT = 0" ■ ■ TTi 

_/l2345\ _/l2345\ _/l2345\ 

^^~1^1 3 2 4 5 )' ^^"1^1 2 3 5 A ) ' ^"yl 5 3 4 2 ) 

We now rewrite the denominator of the right-hand side of equation (14) as 

A{xi, ...,xn)^ (-)''" A(x^(i), . . . , x^N)) = (15) 

A(x^(i), . . . , X^(M)) A(x^(M+l), ■ ■ ■ , X^N)) n i^-^ii) ~ ^<3)) 

ie[i,...,M], je[M+i,...,Ar] 

and observe that 

Ai-l+iV \m-M+N 
—y '-^ ^ {Xi...Xm) Sx{Xi,. . . ,Xm) (16) 

c; A(a;^j(i), . . . ,a;7ri(M)j 

TTifc Dm 

and 

M-l+N 

^ -^MM+l) ■■■•^n2{N) ^ ^ ^^^^ 

C A{Xt,^(M+1), ■ ■ ■ ,X-„2{N)) 

7r2G Jn-m 

After the permutations of type tti and 112 have been performed in equation (14) we remain 
with the following sum over permutations a G Sjv/Sm x Sat-m^ 

Sx{x„...,Xn)^ Y: • • • ^^(M?) ^a(x.(1) ■ ■ ■ X^^M)) _ ^^g^ 

^c: /C n (^CT(i) ~ ^(T(j)) 

o-e Ojv/JMXOiv-M ie[i,...,M], je[M+i,...,Ari 



r(l) ■ ■ ■ •*'7r(M) •^7r(M+l) " " " •*'7r(Ar) 
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E 

<t6 Sjv/SmxSj\ 



AT-M \ 

n^x^J) j SA(a;(^(i) . . .x^(Af))A(a;^(i) . . .x<^(M))A(a;(^(M+i) • • •a;(^(Ar)) 



A(a;(^(i) . ..Xa{N)) 



We insert expression (18) to the sum in equation (12) and apply the Cauchy-Littlewood 
formula. It gives 

^ = (19) 



M N 

n n (1 - XiVj) 
j=i i=i 



E 



I 

M , . 

n (1 - Xa{i)y3) 
,i=i 



A (3:^(1), ■ ■ ■ , '^uiM)) A (a;CT(M+i), ■ ■ ■ , 'Xc(n)) 
A (xo-(i), . . . ,Xo-(Ar)) 



\i,j=l ' / 

It is easy to see that the formula above is equivalent to the formula (11). 

3 Multi-variable integral representation for Cauchy 
transforms of orthogonal polynomials 

In this section we compute (^Z^^[e, H]^ ^, dimif = N, and show that this average is 
equal to jN~ihN-ii^)- Here 7Ar-i is given by equation (10) and /iAr_i(e) is the Cauchy 
transform of the monic polynomial 7r7v-i(ic) (equation (9)). The aim is just to illustrate our 
approach to the computation of the general correlation function )Cn{^, fi) of characteristic 
polynomials on the simple example L = 0, M = 1. 
Once M — 1 algebraic identity (11) takes the form: 

e^"^ j2 ( ^^(1)' \ A(a;^(2),...,x^(jv)) ^20) 



2n[^,H] ,c I e - ^ct(i) / A(X(^(i), . . . ,X(j(jv)) 

Averaging the expression above with the probability measure V'^'^^xi, . . . ,xiy) (equation 
(3)) we notice that each term yields the same contribution {V^^\xi, . . . , xat) is symmetric 
with respect to the permutations). Thus we have 



e 



/ A(x2, ■ ■ ■ , xn)A(xi, . . . , XN)e dxi...dxN (21) 

J e — Xi 



H 

N 



Z]\j J e — X\ 
The integral above can be further rewritten as 

^iV-lg-W(xi) 



dx, ^ dim // = TV - 1 (22) 
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Using the formula (4) we then obtain 

[6, H])^^Ne -—-J ^-^ (23) 



Now we note that 



+ — : (24) 



The first term above is a polynomial of degree A^— 2 which is orthogonal to the polynomial7rAr-i(a;). 
Thus only the second term contributes to the integral. It gives 

{Z],%H])^ = ^M-ihM-i{e) (25) 

or 

AT 

7Ar_i/iAr_i(e) = — / J|(e - Xj) l\'^{xi, . . . ,XN)dxi . . .Xn (26) 

J j=i 

N-l 

where we have used formulas (9), (10) and have expressed Zjy as A^! H c^- 

j=0 

Wc expect that the representation of Cauchy transforms of orthogonal polynomials as 
multi-variable integrals (26) might be known, but we failed to trace it in the standard 
monographs on the subject. 



4 Derivation of the formula for /CAr(e,/i) 

Let us now consider the general case. The correlation function /CAr(e, /i) defined by equa- 
tion (6) is an integral over variables Xi, . . . ,X]\f (which are eigenvalues of H ) with an 
integrand symmetric under permutations. It then follows that each component of the sum 
in equation (11) gives the same contribution to the correlation function ICn{^,P^)- Thus 
we have 

TV' ^ r 
" {N-M)\M\ n J dxi... dxNV^^\xi, ...,xn) 

L 

M ^N-M \ 



X 



U^]tF^ (27) 



1=1 j=M+l 

We decompose the eigenvalue probabihty density function as 

ry ry M N 

'N-M^M 



V^^\x,,...,x^) = ^^^^^P^\{ n (^.-^.f (28) 

i=lj=M+l 

X 7'W(xi, . . . , XM)P(^-^H^M+l, . . . , Xiv) 
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With the decomposition above the integral expression for the correlation function /Cjv(e, /t) 
can be rewritten as 

r rA ( Tl Zn-mZm 



{N -M)\M\ 

j dXi... dXMV^^\xi, . . . ,XM)Fe,[,{xi, . . . ,Xm)I^^~^\xi, . . .,XM,fl') 

where we have introduced two functions of M variables xi, . . . , xm- The first function is 



f,,,{x,,...,xm)= n n ' , (30) 

i J=l fe=i 

The second function, I^^~-'^\xi, . . . , xm, P), can be understood as an averaged value of the 
products of characteristic polynomials over an ensemble oi N — M dimensional Hermitian 
matrices H with eigenvalues xm+i, ■ ■ ■ ,xn, 

7(^-^)(xi,...,XM,/i)= (31) 

/ML _ _ \ 

( n n ^^n-mIxi, H]ZN-M[i^k, H] \ 

\i=l k=l I H 

It is convenient to introduce L + M-dimensional vector f whose components are the 
integration variables xi, . . . , xm and the elements of the vector p,, i.e. 

f = (xi,...,xm,A«i,---,A«l) (32) 

We exploit formula (5) which for l'^^~^\xi, . . . ,xm, P) = l'^^~^\r) gives the following 
expression 

I^''-''\r) = ^ det K-M+,-i(n)],<,,<^+M (33) 
The Vandermonde determinant A(f) can be factorized as 

{L M ^ 
n l[i^j-pik)\^{xi,...,XM)A{fx^,...,fXk) (34) 
fe=ij=i J 

Now we insert equations (33) and (34) to the formula for the correlation function /Cjv(e, fi) 
(equation (29)). It gives: 

where Jm{^, P) is the M- fold integral 

JM(e,/t)= (36) 
jdx,... dx^ det [$.(0:,)],,,,,^ det K_M+.-i(r.)] 



where we have introduced M functions ^i{x) 



M 

n {ej - x) 



(37) 



Thus the computation of the correlation function )Cn(^, P) is reduced to that of the integral 
of M variables. This integral can be further transformed to a determinant form. To 
proceed we first note that det [^^(xj)] can be simplified^ 



det [^i{xj 



A(£) 



(41) 



Insert the expression (41) to the integral (36), rewrite the determinants as sums over 
permutations and find in this way that 



JM(e, A) 



M! 



A(6)A(/i) 



det 



'Kn-M\(^1) 

t^n-m{^m) 



■ T^N+L-l{^l) 

■ t^n+l-i{^m) 

■ 7rAr+L-i(/^i) 



(42) 



where 



^-NV{x)^N-M 



e — X 



T^N-Mil^h) ■ ■ ■ T^N+L-lil^L) 



TTk{x)dx, k&[N -M,N -M + L-l\ 



(43) 



Using the orthogonality of monic polynomials Tiki^) with respect to the measure e ^^'^^^dx 
we observe the relation between -7ffc(e) and the Cauchy transforms hk{e): 



7ffe(e) = -27rie^-^/ife(e) 



(44) 



We insert expressions (42), (44) to formula (35) for the correlation function /Cjv(e, A)- 
Using the relations between coefficients Zj, Cj, 7^- we finally prove (8). 



^This follows from the algebraic relation 



M 



M 

1=1 



= i-r E 



^1 n (efc - eO Xj - Ck 
l^k 



(38) 



which can be derived from identity (11). With the equation above we have 
det [Mxj)] = {-r'det (^) det f e}"^ U 

Noting that 

<-)--(.r'n^) = 4 



(39) 



(40) 



we obtain equation (41) 
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5 Correlation function /C7v(e, jl) and the Riemann-Hilbert 
problem for orthogonal polynomials 

Our formula (8) enables us to express the correlation function /Civ(e!A) i'^ terms of so- 
lutions of Riemann-Hilbert problems for orthogonal polynomials proposed by Fokas, Its 
and Kitaev [15, 16]. It follows from (8) that the correlation function /CAr(e, /i) is deter- 
mined by monic orthogonal polynomials and their Cauchy transforms. In turn, the monic 
orthogonal polynomials and their Cauchy transforms can be understood as elements of 
a (matrix valued) solution of the following Riemann-Hilbert problem. Let contour S be 
the real line oriented from left to right. The upper side of the complex plane with respect 
to the contour will be called the positive one and the down side- the negative one. Once 
an integer n > is fixed the Riemann-Hilbert problem is to find a 2 x 2 matrix valued 
function Y — Y'^'^^{z) satisfying the following conditions: 



Here Y^\z) denotes the hmit of Y^'^\z') as z' ^-^ z E T, from the positive/negative side. 
As is proved by Fokas, Its and Kitaev [15, 16] the solution of the Riemann-Hilbert problem 
is unique and is given by 



On comparing the formulae (8) and (45) we observe that the correlation function ICn{^, fi) 
can be expressed in terms of the matrix elements of the solution of the above Riemann- 
Hilbert problem. The relation provides us with a possibility to investigate the asymptotics 
of /CAr(e, jl) at large for essentially any potential function V{x) entering the probability 
distribution (3). The details will be presented in a forthcoming publication [24]. Here we 
just outline the main steps. First we will show that the formula (8) can be rewritten as a 
determinant whose entries are two kernel functions. Those kernel functions are express- 
ible in terms of the solution of the Riemann-Hilbert problem for orthogonal polynomials. 
The large N asymptotics in the Dyson scaling limit can be studied by the steepest de- 
scent/stationary phase method for Riemann-Hilbert problems introduced by Deift and 
Zhou [19] and developed further in [20, 21, 22, 23]. As a result wc will prove the univer- 
sality of various quantities related to /CAr(e, fi). In particular wc shall be able to prove the 
universality of the distributions of the level curvatures and the local density of states, as 
well as to derive the Poisson kernel distribution of S-matrix in a random matrix model of 
quantum chaotic scattering. Understanding of universality of objects arising in random 
matrix theory in various scaling limits was recently a subject of intensive work both in 
physical [26] as well as mathematical [27, 17] communities. 



• Y^"\z) - analytic in C \ E 






(45) 
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6 "Duality relations" for the Gaussian case 



For the particular case of the Gaussian measure V{x) = equation (5) for the ensemble 
average of products of characteristic polynomials can be written in an equivalent form: 

/ ci[^]e-^'^^^' det (a (8) Iat - 1l <8) ^) oc / [dQ]e- det (X - iQ^ (46) 

Here the integration in the left-hand side goes over the manifold of Hermitian matrices of 
the size N x N, whereas in the right-hand side it goes over L x L Hermitian matrices. 

This integral identity can be viewed as a certain "duality" relation for matrix inte- 
grals. It emerged in various physical contexts, most notably in the context of matrix 
models of the string theory (see e.g. [28], p. 27) where it played an important role in 
understanding equivalence between one-matrix models of the quantum gravity and the 
so-called Kontsevich-type models [28, 29]. To understand equation (46) one recalls that 
the orthogonal polynomials for the Gaussian case are Hermite polynomials possessing an 
integral representation: 

7r,(A) oc e^^V2 T dqq'e-''('^-''') (47) 

J—oo 

Substituting (47) into (5) one easily brings the right-hand side of the latter to the form 
n Zn[X,,H]\ oc ^ / rfg.A{g,}e-^[^^^^^^^^^] detQ^ (48) 



\j=l I H 



where Qx = diag(gi, . . . , qi)- From the other hand the expression above can be obtained 
after shifting X — iQ ^ —iQ in equation (46), diagonalizing the L x L matrix Q and 
integrating out the angular degrees of freedom with help of the Itzykson-Zuber-Harish- 
Chandra formula. A particular case of the formula (48) was mentioned in [5] in the context 
of symmetric Jack polynomials. 

It is easy to check that the Cauchy transforms of the Hermite polynomials are given 

by 

h,{e) <x J dqq'eA"^-''^^^'^'^'^) (49) 

Here the integration domain is < g < oo for Im e > and — oo < g < for Im e < 0. 
This gives us a possibility to rewrite our main object - correlation function (8) - in the 
form of M -|- L-fold integral [30] : 

ICNii, A) oc ^lly^^^y / dQA{Q}e-''[^^^'-'^^'^] det g^"^ (50) 

where Q — {qi, ■ ■ ■,qL, Ql+i, • • • , Qm+l) and E — {jl, e+, — e_), with e+ and e_ denoting 
spectral parameters efc with positive (negative) imaginary part, respectively. The inte- 
gration domain is — oo < < oo for 1 < Z < L but < < oo or — oo < g^ < for 
L < I < M + L, depending on the sign of the imaginary part of the corresponding spectral 
parameter. 
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Expression (50) is a generalization of duality relation (48). It arises most naturally in 
the method based on Gaussian integral representations and Itzykson-Zuber type integra- 
tions [11]. Similar identities hold for the chiral ensemble (Laguerre polynomials), see [13] 
for more details. 

7 Conclusions 

In this paper we have found an exact formula for the general correlation function ]Cn{^, ft) 
containing both products and ratios of characteristic polynomials. Our result is valid 
for an arbitrary ensemble of Hermitian matrices of (3 = 2 class. The obtained formula 
establishes a correspondence between the correlation function /Cjv(e, A) and the Riemann- 
Hilbert problem for orthogonal polynomials. It is remarkable as it enables us to study 
the large N asymptotics of this correlation function via the Riemann-Hilbert approach. 
Among interesting prospects for future research we woTild like to mention a challenging 
problem of extending our calculations to other symmetry classes (3 = 1,4 (cf. [31]) as 
well as to ensembles of non-Hermitian random matrices important for the problems of 
quantum chaotic scattering [32]. 
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